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Abstract. This paper discusses new perspectives and approaches to the problem of disk dynamics where, in this 
study, we focus on the effects of viscous instabilities influenced by boundary effects. The Boussinesq approximation 
of the viscous large shearing box equations is analyzed in which the azimuthal length scale of the disturbance is 
much larger than the radial and vertical scales. We examine the stability of a non-axisymmetric potential vorticity 
mode, i.e. a PV-anomaly. in a configuration in which buoyant convection and the strato-rotational instability do 
not to operate. We consider a series of boundary conditions which show the PV-anomaly to be unstable both on a 
finite and semi-infinite radial domains. We find these conditions leading to an instability which is the disk analog 
of Tollmien-Schlichting waves. When the viscosity is weak, evidence of the instability is most pronounced by the 
emergence of a vortex sheet at the critical layer located away from the boundary where the instability is generated. 
For some boundary conditions a necessary criterion for the onset of instability for vertical wavelengths that are a 
sizable fraction of the layer's thickness and when the viscosity is small is that the appropriate Froude number of the 
fiow be greater than one. This instability persists if more realistic boundary conditions are applied, although the 
criterion on the Froude number is more complicated. The unstable waves studied here share qualitative features 
to the instability seen in rotating Blasius boundary layers. The implications of these results are discussed. An 
overall new strategy for exploring and interpreting disk instability mechanisms is also suggested. 
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I> 1. Introduction Keplerian profile (Li et al. 2000), baroclinic instabilities 

(Klahr & Bodenheimer, 2003, Johnson & Gammie, 2005 

The maffneto-rotational mstability (Balbus, 2003, MRI ^ , i i a o onn-7\ j. • j. j-i, j j. • j 

J, , ° ,,,,,, Petersen et al. A-B, 2007), transient growth and sustamed 

!3 ■ heremaiter) IS generally considered to be the leadme can- u -i- i j • /d- u j o r, i mnn t p 

, , . . , r , , , subcritical dynamics (Richard & Zahn, 1999, lounnou & 

■ - ■ didate explaining the source of enhanced transport ob- ^akouris, 2000, Chagchshvih et al. 2003, Tevzadze et al. 

served for disk systems. Three conditions are required 2003, Yecko 2004, Umurhan & Regev, 2004, Umurhan et 

for its operation: the concurrent presence of rotation and , onoo td p -mi onnc t p t j-j.- 

^ , „ , , 3.1. 2006, Barranco & Marcus, 2005, Lesur & Longaretti, 

shear, a primordial (no matter how small) magnetic held, onoc t -iu -1 000-7 j ii, • i ^d- j 

' . . . \ n 1 n -1 1 1 • • 1 2005, Lithwick, 2007, and see the experiments 01 Richard, 

and sufficient ionization 01 the nuid so that the gas IS m the r,nni j t- j. 1 nnr,/i\ j j j 

. , . . , T 1 2001 and Ji et al., 2006) and, unsteady wave dynam- 

MHD regime. Consequently, it IS natural to ask the ques- ■ ti n 1 ■ n ■ 1 ■ ^ u-i-j. r-r, 1 

, . . ,. , ^ ^ ICS like the Fapalouizou-Prmgle instability (Papaiouizou 

tion: What happens man accretion disk if one or more of ^ p^.^^^^^ ^g^^^^ hereinafter "PPI", and the strato- 

these conditions are not satisfied? , ,. 1 • . 1 -i-i /t-> u n j. 1 000c ci, 1 ui p 

rotational instability (Dubrulie et al., 2005, Shalybkov & 

Supposing that there are disks in which the MRI or any ^^diger, 2005, Umurhan, 2006, Brandenburg & Riidiger, 

other MHD (dynamo) mechamsm is either weakly operat- 20O6) , hereinafter "SRI" . The majority of these recent in- 

ing or entirely absent: what else can drive activity, possi- ^egtigations (excepting Yecko, 2004, and Afshordi et al. 

bly even leading to turbulence? Attempts to answer this 2005) have focused on strictly inviscid processes. Indeed, 

question include, but are not limited to, defects in the j.i i • 1 i j. r, • j. • j.- j. 

^ ' ' the classical approach to such questions is to investigate 

Send offpnnt requests to: O.M. Umurhan, e-mafl: processes which might lead to turbulent transport by first 

o.uinurhan@qmul.ac.uk establishing a mechanism of hnear instability from the 
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Since astrophysical fluids have some effective viscosity 
- however small it may be - we pose the question: could a 
weakly viscous flow in a sheared and rotating environment 
undergo an intrinsically viscous type of linear instability 
that nonlinearly saturates with significant amplitude? A 
few previous studies have addressed this question. Kato 
(1978) demonstrated that if a fluid's viscosity is a func- 
tion of the state of the fluid then a disk can experience 
pulsational dynamics in a way similar in quality to stel- 
lar pulsation like that in the theory of Cepheid variables. 
Hereafter we shall refer to this effect as the viscous pulsa- 
tional instability (VPI). Latter & Ogilvie (2006) reexam- 
ined the VPI by studying how even axisymmetric f-modes, 
in a shearing sheet environment, create fluctuating stresses 
that explicitly draw energy from the shear which leads to 
overstability. Kleiber & Glatzel (1999) have shown that 
accretion tori (in this case, ones which have a constant 
specific angular momentum distribution) can be linearly 
unstable above a minimum Reynolds number. DubruUe et 
al. (2005) and Shalybkov & Riidiger (2005) also report 
that the growth rate of the SRI may be, under certain 
conditions, enhanced by viscosity. 

Before proceeding we present some remarks concern- 
ing the nature of the fluid state to bo studied as well as 
a discussion about boundary conditions and a proposed 
alternative way to consider their uses. 

1.1. Viscous considerations 

The non-commutative nature of the Navier-Stokes equa- 
tions in the limit of Rc^ oo with the Eulcr equations is 
a long standing fact (see the discussion in Schlichting & 
Gersten 2001, pg 96-8). For instance, viscous stresses do 
not necessarily vanish on the boundaries of a viscous flow 
when the viscosity limits to zero. In turn this implies that 
one may not properly take the Navier Stokes equations 
and naively substitute zero viscosity to reach the inviscid 
limit (Lions 1993, Joseph 2003). Aside from the generic 
appearance of boundary layers, other effects can appear 
when a weak viscosity is included into problems of study. 
For example, in the nearly inviscid shallow-water theory of 
strong shear flows (Balmforth, 1999), normal- modes can 
emerge out of a continuous spectrum when viscosity is 
introduced into the dynamics. 

Another example of the subtleties inherent in vis- 
cous flow is the instability associated ToUmien-Schlichting 
waves (see discussion of T-S waves in Schlichting & 
Gersten, 2001, and Schmid & Henningson, 2000) - which 
are traveling waves appearing in wall-bounded flows that 
neither grow nor decay in the inviscid limit and become 
unstable when viscosity is included in the analysis (ex- 
amples include plane Pouiselle-flow and Blasius bound- 
ary layers). Similar to the process discussed by Balmforth 
(1999), the traveling mode becomes unstable through the 
interaction of an inviscid normal mode and a viscous nor- 
mal mode - the latter of which exists only as a member 
of the continuous spectrum in the inviscid limit (Baines, 



Majumdar & Mitsudera, 1996). Far from being considered 
mathematical oddities, T-S waves appear to play a promi- 
nent role in the transition to turbulence in boundary-layer 
flows (for a recent summary see Drazin, 2002). A situation 
studied which closely resembles the condition of an astro- 
physical disk is the modal and non-modal response of a 
rotating Blasius boundary layer (Yecko & Rossi 2004) in 
which instability is promoted when the azimuthal scale of 
a disturbance is longer than its vertical scale. 

In a general sense, because the governing equations are 
of a higher order in the viscous case, a new space of pos- 
sible solutions emerges which are either absent or inactive 
in the inviscid case. Our main query is therefore: if rota- 
tionally supported flows are (linearly) well-behaved in the 
exactly inviscid limit but the viscous flow shows some type 
of dynamically signiflcant behaviour - even as the inviscid 
limit is approached - then might it be misleading to test 
stability only of exactly (Re = oo) inviscid flows? Perhaps 
the subtle nature of disks is linked to this feature. That 
which best summarizes this perspective is the quote at- 
tributed to the atmospheric dynamicist E. T. Eady where 
he is purported to have said, "It is not the process of lin- 
earization that limits insight. It is the nature of the state 
we choose to linearize about," (Bayley, Orszag & Herbert, 
1988). 

To prospect for an instability mechanism that might 
lead to sustained unsteady behavior by assuming a turbu- 
lent viscosity model a priori might seem contradictory at 
first. But, given the difference in behavior known to exist 
(in other problems) between nearly inviscid and exactly 
inviscid models it is therefore mandatory to clarify these 
differences within the context of an astrophysical disk too. 
A conjecture that such an investigation could address is 
the following. It is reasonable to suppose that disks are 
continuously fed with a turbulent flow fleld either by in- 
fall or some mass transfer processes. Could this turbulent 
flow field undergo a secondary transition into another dy- 
namical state (possibly turbulence of a different stripe) 
due to the turbulently enhanced viscosity? Suggestions 
which hark on these lines of thought are found in Doering 
et al. (2000), Kersale et al. (2004) and more directly in 
terms of scicondary transitions induced by Ekmann flow 
as suggested in Lesur & Longaretti (2005). Overstability 
driven by material fluctuations in the turbulent stresses 
(i.e. the VPI, Kato, 1978 and Latter & Ogilvie, 2006) are 
also candidates for such secondary transitions. 

1.2. Interpreting boundary conditions and their effects 
on dynamics 

Much attention has been devoted to evaluating processes 
which are intrinisic in some way to the fluid - mean- 
ing to say that it is assumed that it is more valuable 
to study those mechanical processes which are minimally 
sensitive to the boundary conditions imposed and max- 
imally "emerging" out of something essential about the 
fluid and its basic state. The MRI is an example of this as 
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well as other more basic fluid dynamical instabilities such 
as the Rayleigh- Taylor and Rayleigh-Benard instabilities. 
We wonder if this approach to the question of disk turbu- 
lence may be self-limiting given that many fluid instabil- 
ities which lead to some form of turbulence in terrestrial 
flows are driven in large part by the boundary conditions 
of the system (e.g. T-S waves and turbulent transition). 
An alternative way is to view boundary conditions as a 
filter for certain solutions or as a tool to classify solu- 
tions. Kersale et al. (2004) adopt a similar philosophy by 
studying the linear response of an incompressible fluid in 
a Taylor-Couette type of cylindrical flow subject to a vari- 
ety of boundary conditions. Of course, the Taylor-Couette 
setup and the boundary conditions they test are not what 
one would expect in a terrestrial apparatus or experiment, 
however, if taken as a metaphor for a disk environment 
then this sort of exploration allows one to test, evaluate, 
understand and eventually categorize the dynamical re- 
sponse of a fluid as a function of differing boundary con- 
ditions. 

The inviscid PPI and SRI are examples of linear in- 
stabilities which come about due to the imposition of ar- 
guably artificial boundary conditions on inner and outer 
walls of a model disk system. On the other hand, these 
results may be intepreted in terms of the Hayashi- Young 
criterion for wave instability which states that a physi- 
cally separated wave pair may promote linear instability 
if the waves counterpropagate with respect to each other 
with nearly the same wavespeed and if the waves have 
an action-at-distance effect upon each other (Hayash i & 
Young, 1987 and see also Sakai, 1989, Baines & Mitsudera, 
1994). Indeed Goldreich, Goodman and Narayan (1985)0 
point out that the PPI may be viewed as a process result- 
ing from the interaction of a pair of edgewaves mutually in- 
teracting with each other across a wave-evanescent region 
H. The SRI may also be similarly rationalized (Umurhan, 
2008). 

Thus although the counterpropagating edgewaves re- 
sponsible for the PPI and SRI are understood to result 
from the use of unrealistic boundary conditions, it is cer- 
tainly not ruled out that the general counterpropagating 
wave mechanism could be at work in real disks. The linear 
instability of disks with two or more (potential) vorticity 
defects (e.g. Li et al., 2000) could be interpreted as an 
instance of this process. 



^ By this it is meant to say that there exists a wave evanes- 
cent region separating the waves in question. 

They show that the interior of the domain does not intrin- 
sically support propagating waves as it behaves like a wave 
evanescent zone. However the imposition of the boundary con- 
ditions brings into existence waves that propagate along the 
radial boundaries of the domain. 

^ Examples of this process are well known in atmospheric 
flow (Charney & Stern, 1965, Hoskins et al., 1985, Davies & 
Bishop, 1994). 



1.3. An overview of the findings in this study 

From a systematic asymptotic scaling analysis we derive 
in Section 2 and Appendix A the equations appropriate 
to a box section of a viscous shearing accretion disk (as- 
suming an a viscosity formalism) by exploiting the small- 
ness of the parameter e which assesses the ratio of the 
soundspeed to the rotation speed measured at some ra- 
dial point of a circumstellar disk. We refer to this model as 
the Viscous Large Shearing Box (VLSB) and these equa- 
tions have appeared before (cf. Latter & Ogilvie, 2006). 
We are reminded that while the velocity fluctuations in 
the shearing box are an order e smaller than prevailing 
rotational ( "Keplerian" ) velocities, the steady accretion 
velocities implied by the alpha viscosity model are an or- 
der smaller than the same disk rotational velocity. 

We consider the fate of a non-axisymmetric potential 
vorticity disturbance (or simply "P V-anomaly" ) subject 
to varying boundary conditions. Accordingly, in Section 
3 and 4 the VLSB are analyzed in the limit where the 
perturbation's azimuthal length scale is asymptotically 
larger than its radial and vertical scales (i.e. the quasi- 
hydrostatic semigeostrophic hmit, Umurhan, 2006, QHSG 
for short). Additionally we assume that the vertical com- 
ponent of gravity and entropy gradient are constants. 

In Section 4.1 we formulate energy integrals of this 
reduced system in order to better understand what can 
contribute to destabilizing the PV-anomaly . The energy 
budget is characterized by a Reynolds-Orr type of equa- 
tion whose sources and sinks are given by the energy which 
the PV-anomaly can extract from the shear, receive from 
the boundaries or lose due to dissipation. 

In the spirit of Kersale et al. (2004) we analyze the re- 
sponse of the PV-anomaly subject to a controlled array of 
boundary conditions. In discussing boundary conditions 
we refer to the boundary closest to central object as star- 
side as opposed to the side furthest away from the object 
to which refer to as farside. We consider the dynamics as 
occurring on both a semi-infinite domain (farside at in- 
finity) and on a finite domain. Below we summarize the 
main findings. Note that we have made sure to consider 
boundary conditions which filter out the SRI or PPI in- 
stabilities. 

In Section 4.2 an asymptotic analysis is done for the 
limit where the scaled turbulent viscosity parameter (de- 
flncd in the text as e) is small. We flnd instability if the 
Froude number of the flow exceeds 1 for modestly large 
vertical wavenumber. Additionally, the PV-anomoly inter- 
acts with a critical layer of the flow creating a potential 
vorticity sheet sheet whose radial extent is e^^^ the size of 
the vertical extent of the disk. This analysis illustrates how 
an inviscid edgewave phenomenon (due to the no-normal 
flow starside boundary condition) becomes unstable when 
viscosity is included. Most importantly is that the insta- 
bility is driven in part by the injection of energy through 
the boundary. 

We consider in Section 4.3 finite domain disturbances 
of the PV-anomaly and let the viscosity parameter be an 
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order 1 quantity. The fourth order normal mode problem 
requires us to appeal to numerical computational meth- 
ods for solutions. We impose on the farside boundary that 
both the disturbance pressures and PV-anomalies vanish. 
At the starside boundary we require that there be no- 
normal flow there (as above). The remaining starside con- 
dition takes on four possibilities: (a) the flow is rigidly 
coupled at the wall, (b) the perturbations are stress-free, 
(c) the PV-anomaly is zero, (d) the PV-anomaly gradi- 
ent is zero. The first two of these conditions are physi- 
cally realistic. The latter two offer a means to consider 
the the effect of energy injection (or lack thereof) through 
the boundaries and to compare with the analytical anal- 
ysis. For rigid and stress-free boundary conditions we see 
clear indications of a ToUmien-Schlichting type of instabil- 
ity, similar to the instability of rotating Blasius boundary 
layers (Yecko & Rossi, 2004) and the energy budget of the 
disturbances show that this process does not draw upon 
energy across the boundaries. 

2. Viscous Large Shearing Box and its QHSG 
approximation 

In Appendix A we consider a box section of an a-disk 
centered about its midplane and at a distance from 
the central object. If the disk is cold, then it means that 
the quantity defined by the ratio of the typical value of the 
local midplane disk soundspeed, Cs, to the local Keplerian 
velocity, Vq, 

Cs 

is less than 1 by some substantial amount: protoplane- 
tary disks, for example, are believed to have an e « 1/20. 
Using now familiar scaling arguments and exploiting the 
smallness of s we derive from the full equations of mo- 
tion in cylindrical coordinates appropriate equations of 
motion in what we refer to as the Viscous Large-Shearing 
Box (VLSB for short). The tactics and procedures behind 
this effort are the same ones employed in the derivation of 
the Large-Shearing Box (LSB) (Umurhan & Regev, 2004) 
however, the viscous stresses are included. We have. 



{dt - qD.oxdy)p -f V ■ (pfc -I- p)W = 0, 
{dt - qVlQxdy)u' + u' ■ Vm' - 2fiou' = 



dxP 



Pb + P 

{dt - qVlQxdy)v' + u' • Vu' + (2 - q)VlQu' 

9yP 
Pb + P 
dzP + pg{z 



(1) 



+ K (2) 



Pb + P 



+ K' (3) 

+ K, (4) 



{dt — qfloxdy)w + u • Vw = - 
{dt - qnoxdy)J: + u' ■ VT, = (5) 
in which the total entropy is defined by 
Pb+P 



and 7 is the ratio of the specific heats at constant pres- 
sure to the specific heat at constant volume. The vertical 
component of gravity is dependent on z 



(6) 



All primed quantities are perturbations about the basic 
flow. The viscous stresses are 

{pb + p)K = V{dl + dy)n' + d^fid^u' 

+dx [r~l{dxu' + dyv') + d^fiw'] - f a,77V • u', (7) 
{pb + p)% - v{dl + dl)v' - q^od.fi + d,fid,v' 

+dy [fi{dxu' + dyv') + d.fjw'] - fayryV • u', (8) 
{pb + p)K - V{dl + dy)w' + d^fid^w' + 

+f,d,{dxu' + dyv') + d,f,d,w' - Id^fiW ■ u', (9) 



where 



P)- 



(10) 



The above equations are non-dimensional. Time is scaled 
by the local rotation time of the box. All lengths are 
scaled according to a length H <^ which is compa- 
rable to the disk thickness (see Appendix A). Pressures 
are scaled according to the product of the local mid- 
plane soundspeed and density, which is in turn based on 
some fiducial characteristic temperature scale. For fur- 
ther details see Umurhan & Regev (2004). x represents 
the radial (shearwise) coordinate of the SB while y is 
the azimuthal (streamwise) and z is the vertical coordi- 
nate (normal to the disk midplane). The velocity compo- 
nents, i.e. u' = {u' ,v' ,w'}, are for the radial, azimuthal 
and vertical directions. It is important to keep in mind 
that these fiow variables represent perturbations about 
the steady Keplerian fiow. Oq, sometimes also referred to 
as the Coriolis parameter, is 1 in these nondimensional- 
ized units, meaning to say because time has been scaled 
according to the dimensional value of the rotation rate 
at Rq, i.e. ri(i?o), the local Coriolis parameter formally is 
equal to one. We retain this symbol in order to flag the 
Coriolis effects in this calculation. The local shear gradient 
is defined to be 



(11) 



in which fl{R) is the full disk rotation rate. For Keplerian 
disks the value of q is 3/2. The local Keplerian fiow is rep- 
resented here by a linear shear in the azimuthal direction, 
i.e. —qfloxy. 

The steady state quantities are denoted with index b 
and in particular we assume that pressure (p^ ) and density 
{p^ ) profiles satisfy the hydrostatic balance relationship 



-R 









9zP, = -P,9{z). 



(12) 



E = In- 



{pb + pV 



All corrections to this equations are of a higher order 
and ignored here. The expression responsible for the VPI 
can be identified as fluctuating viscosity parameter in ([8|) , 
term —qilodxij- 



5 



We call to attention that the accretion and meridional 
velocities characterizing a-disks (Kluzniak & Kita, 1999) 
do not appear in the VLSB equation set (see Appendix A) . 
By definition, turbulent disks exhibit accretion velocities 
as they are the natural consequence of equations describ- 
ing global dynamics. The shear velocities are in fact quite 
complex as Kluzniak & Kita (1999) showed for the par- 
ticular case of an Shakura-Sunyaev type of a-disk. The 
radial velocities in steady state are found to be sheared 
in the vertical direction and, as well, there exists a verti- 
cal component to the flow with both radial and vertical 
dependence (the meridional flow) . However the scaling ar- 
guments implemented to reach these "shearing box" equa- 
tions, especially the relative scaling relationships between 
the dynamical velocities and the accretion scalings, show 
that the influence of the steady accretion rate appears at 
higher orders in the expansion procedure. In other words, 
dynamical perturbations on the scale of the box do not feel 
the effects of steady accretion and meridional flow - they 
only feel the effects of the steady Keplerian shear. The 
scaling analysis also shows that the a-viscosity (which is 
the driver of the accretion flow) does influence the dynam- 
ics at these scales and is the reason why it appears in these 
equations. 

3. The quasi-hydrostatic semigeostrophic 
approximation of the VLSB equations 

The equations of motion may be simplified for fur- 
ther analysis by implementing the quasi-hydrostatic semi- 
geostrophic (QHSG) scaling arguments used in Umurhan 
(2006). The QHSG is useful in its ability to expose 
the essential mathematical features of the inviscid-SRI 
(DubruUe et al., 2005, Umurhan, 2006). 

We suppose that the azimuthal scales of motion are 
much larger than the radial or vertical scales. We mea- 
sure this with the small parameter 5. In order to maintain 
asymptotic validity we assume the following orderings 



£ < (5 < 1. 



(13) 



Thus we suppose that the following operations upon 
dynamical quantities scale accordingly as 

Then we suppose that the radial and vertical velocities are 
correspondingly smaller than the azimuthal velocities by 
this same scale, in other words 

v'-O(l), u\w'^0{5). 

These scalings say then that 

—qx^L^dy -\- u'dx + v'dy + w'dz ^ O {S) . 

These scalings will make it easy to follow waves propagat- 
ing with respect to the background Keplerian flow velocity. 
Therefore, the temporal dependence should also scale by 
the scaling appropriate to dy. It follows that 



Furthermore we say that the density, pressure (and by 
implication, the entropy) fluctuations are all order 1, that 
is 

p,p,S~0(l). 

The new issue that must be addressed here is to suggest 
a scaling that brings in the viscous terms at the lowest 
non-trivial order. To this end setting a ^ O {S) achieves 
this goal and we shall formally write a — dai. In sum, 
then, to lowest order we have the following reduced set: 

{dt-qnoxdy)p + \7 -{pb+pW ^0, (14) 
dxP 
Pb + 

{dt - qnoxdy)v' + u' • Vi/ + (2 - q)nou = 

Pb + P ^ ' 

= -dzp - pg{z), (17) 
{dt-qnoxdy)Y. + u! ■SJiY.b + Y.)^^, (18) 

where we have introduced the basic state entropy St, and 
its dynamically varying counterpart S which are defined 

by 

Vb 1 + - 

^, S = ln^ 2^7- (19) 

Pb 



= 2VLnv' 



(15) 



T,b = In 





\ T 




--) 




Pb ) 



Only the azimuthal direction stress component survives at 
lowest order due to this scaling argument, 

{pb + p)N'y = f\dlv' - qQodxV + dzfjdzv'. (20) 

Although we have invoked scaling argments leading to the 
above sets of equations we have not formally rewritten all 
of the variables to signify these assumptions as it is our 
desire to preserve the transparency of the subsequent pre- 
sentation. Note that effect responsible for the VPI sur- 
vives this scaling argument as it appears in (j20|) as the 
term —qflodxfi- 

4. Boussinesq Simplification, Assumptions and 
Linearized Dynamics 

In Umurhan (2006) it was shown that the QHSG ap- 
proximation of Boussinesq disk models recovers the lin- 
earized hydrodynamic behavior contained therein for con- 
current small values of the azimuthal wavenumber and 
wavespeed. It was further demonstrated that the dynam- 
ics contained in the QHSG approximation of the LSB is 
faithfully represented if one considers instead the equiva- 
lent incompressible Boussinesq (Spiegel & Veronis, 1960) 
version of QHSG approximated LSB equations. Applying 
this sequenced reasoning to the linearized version of (|14t - 
fT51) gives, 

dxU + dyV + dzW = 0, (21) 
= 2w~dxTi, (22) 
{dt - qxdx) v+{2-q)u ^ -dyli + N'y, (23) 
= -dzW + Q, (24) 
{dt-qxdx)Q = -N^w, (25) 
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where we have dropped all primes from the velocity quan- 
tities and explicitly set fig to its value of 1. In the usual 
Boussinesq approximation, density fluctuations are dy- 
namically significant when coupled to gravity. In these cir- 
cumstances p is replaced hy —6. The non-dimensionalized 
temperature quantity 9 and its associated steady state 
temperature field Tb{z) are characterized by the (lin- 
earized) conservation relation, 

{dt - qxdy)e + wd.Tb = 0. 

For clarity we reexpress this thermal quantity in terms of 
Q given hy Q = gd/ pb^. In the Boussinesq approximation 
fluctuating density variables influence the dynamics when 
coupled to gravity. Forthwith, pb and Pb are taken to be 
constant (set to 1) and, furthermore, pb is absorbed into 
the fluctuating pressure leading to defining the enthalpy 
n = p/pb- AH unprimed velocities (i.e. u,v,w) are now 
understood to represent linearized disturbances. The non- 
dimensionalized Brunt- VdisaJia frequency^ N, emerges in 
the equation for the perturbation temperature field (|25p 
and is given by 



N 



2 _ 



gj;d.Tb. 



(26) 



Throughout this study N is taken to be real (buoyantly 
stable). The azimuthal stress is 

N'y^^[{dl+d'z)v~qdxU\. (27) 

The term that gives rise to the VPI appears in the above 
as ~ —qdx^- We proceed further by (i) operating on (p3)) 
with dx, (ii) operating on (j22p with dy, (iii) and subtract- 
ing the results to reveal 

{dt - qxdy)dxv = (2 - q)d,w + dxN'y, (28) 

where the incompressibility condition was used in writing 
the first term on the RHS of this expression. Multiplying 
(PSj) by (2 — q)/N'^ followed by operating on the result 
with dz gives 



{dt - qxdy){2 ~ q)d,{Q/N^) = -(2 - q)d. 



(29) 



Adding these two equations together and making use of 
the relationships (|22p and yields the following single 
equation for 11: 

{dt - qxdy) {dl + d^F^d,) n = 

l-{dl + d^^)dlu-l-2qdln. (30) 

To showcase the elements in the above we have written 
2ai/3 ■s-^ Re~^ in order to remind ourselves that the ai 
parameter is like the inverse of a Reynolds Number. We 
have also defined the epicyclic Froude number 



f: 



_ 2(2 -g) 

- 



* A momentary comparison to (|f 8|) should convince the 
reader that represents a perturbed entropy quantity 
(DubruUe et al., 2005) making (|25|) a reasonable analog of the 
linearized form of ([5]). 

^ However, note that the actual Reynolds number, denoted 
by R, is a factor of 1/S larger than Re, i.e. R = Ke/S. 



which is in general a function of the vertical coordinate 
z and vanishes on the symmetry axis driving the local 
Froude number to very large values. As we shall see, F^ 
determines the character of the solutions that emerge. 
Umurhan (2006) demonstrates by a direct comparison of 
calculations that the onset of the SRI (in the QHSG limit 
of the LSB equations explored there) is reasonably well 
captured analytically when TV^ and g are assumed to be 
constants. Guided by these previous results as well as sim- 
ilar use in a series of other studies (e.g. Tevzadze, et al., 
2004, Bodo et al. A-B, 2007) we shall assume 



g — constant, N — constant, 



(31) 



implying that i^c is a constant as well, and we restrict our 
considerations to vertically periodic solutions. We shall 
return to this matter in the Discussion. 
In this asymptotic theory 

Q = {dl + d^Fld,)^ 

corresponds to a the perturbed potential vorticity 
(Tevzadze et al. 2004, Umurhan, 2006), also known as the 
potential vorticity anomaly (Hoskins et al., 1985) and also 
called vortensity in the astrophysical literature (e.g. Klahr 
& Bodenheimer, 2003). Throughout the rest of this work 
we will interchangeably use the terms potential vorticity 
perturbation, potential vorticity disturbance and potential 
vorticity anomaly (i.e. PV-anomaly). The PV-anomaly Q 
relates to a vorticity pointing in the vertical direction. 

We consider travelling wave normal mode solutions to 
(|30p : the modes are assumed to be azimuthally periodic 
on scale Ly and vertically periodic on scale Lz, 



n = n(a;)e''=('?'=*+«)+*'^^ + 



c.c. 



(32) 



where the azimuthal wavevector k can be any positive 
number while the vertical wavevector /3 is any real number. 
The wavespeed c can be complex: when Im(c) < there 
is growth of the wave. Furthermore we define 



1 



2ai 



3kq Kekq 



(33) 



to be a wavenumber scaling of the viscous parameter. Thus 
the equation governing the structure function Il{x) be- 
comes 



i{c - x){dl -kl)tl = e{dl -p^ - 2q)dltl. 



(34) 



The Froude-wavenumber is defined as k^ = F^p. ([M]) will 
be the fundamental equation of study. Furthermore ex- 
pressed in this form, the azimuthal and radial velocities 
are 



^dxTl, 



2-q ^ 
qk 



i(c-x)z9,n-i(i)n 

+Hdl -p^- 2q)dxX 



(35) 



(36) 
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with the vertical velocity following from evaluating the in- 
compressibility equation at the boundary and using ([34]), 



I.e. 



w 



qkkj, _^ 



2/3(2 -qy 



(37) 



Finally, the perturbed potential vorticity is Q = {d^ — 
kp)Il. The remainder of the boundary conditions, namely 
in the x direction, will be stated in the following sections 
according to the problem being solved. 

4.1. Energy Integrals 

It is instructive to develop global energy integrals as such 
quantities aid in developing an interpretation of the re- 
sults in the following sections. We restrict our attention 
to the energetics associated with the QHSG-Boussinesq 
model set (I22j|2ip keeping in mind the assumptions we 
made about vertical and azimuthal periodicity and the 
constancy of iV^. The radial conditions are left arbitrary 
and they will be dealt with accordingly in each of the 
subsequent subsections. We proceed by defining the per- 
turbation thermomechanical energy density £ as 



£ = 



1 



iV2 



by multiplying (gS]) by w and by Q/N^ and inte- 
grating over a domain which is periodic in the vertical 
and horizontal directions, and finite in the radial direc- 
tion given by < a: < a:i, where xi is left arbitrary, we 
find 



dE 
'dt 



-^shcar ~t~ £^strcss -^visc; 



(38) 



in which volume integrals are 



E= / £d^'x, 



V 



E'shcar = 9 / X, 
{\dM^ + \dM^ + 2q\v\^)d^^. 



(39) 



integrated on volume V with the volume element d'^x 
dxdydz and where the surface term is 



7] 

Hu H vdxV 

Pb 



dydz. 



(40) 



In writing iSgtress we have made use of the periodicity 
conditions in the vertical and horizontal directions. The 
bracketed term as appearing means 



/(Xout) - /(O), 



where Xout is the location of the outer boundary. 

The energy integral ([55)1 is the Reynolds-Orr Equation 
appropriate for this QHSG system. The energy E is com- 
posed of the baroclinic thermal term (~ 9^) plus a kinetic 



energy term (^ 



however, the kinetic energy term con- 



vertical and horizontal velocity contributions are small by 
comparison in the scaled system of equations according to 
the QHSG approximation of the original set (see Section 
3 and Umurhan, 2006, for details). The term —Evisc rep- 
resents the integrated losses and is comprised of viscous 
losses due to the azimuthal velocity since, by the same 
reasoning as above, the corresponding losses due to verti- 
cal and radial velocities (in the scaled system) are negli- 
gible. The perturbed viscous stress term, responsible for 
the VPI plays, a destabilizing role for these PV-anomalies 
as it appears in the above as the the term proportional 
to 2g|wp. The total external stresses on the system is 
given by i^stress and is comprised of the surface integrated 
body pressure and the surface viscous stress - the latter of 
which is expressed only with the tangential stress due to 
w. Finally E'shear is the Reynolds stress due to the back- 
ground shear state. This expression may be interpreted as 
accounting for the amount of energy perturbations extract 
from the background shear state. We note also that in- 
spection shows that i?visc > always, while the remaining 
terms i?stress and -Eshcar may be either positive or nega- 
tive given the state of the perturbed flow or the boundary 
conditions employed. 

4.2. Asymptotic theory on a semi-infinite domain 

The following assumptions are made in order to proceed 
analytically: (i) the domain in the x direction lies between 
and oo, thus we have Xout ~ oo, (ii) we require that all 
quantities decay as x ^ cxd, (iii) there is no normal- flow 
at X = - this inner location is considered starside as it 
is closest to the central object, (iv) the viscosity is weak 
but finite, hence, we assume that e ^ l[£l (v) even though 
there is viscosity in the problem we impose no particular 
stress condition at x = and let the fluid quantities be 
dictated by what emerges in the interior of the domain. 
This essentially means that flow stresses at the starside 
are allowed to adjust according to the dynamical response 
happening in the interior of the domain. In practice it 
translates to only enforcing the no-normal flow boundary 
condition^ 

4.2.1. Expansions and outer solution 

We assume the following expansions well aware that this 
is a singular perturbation calculation because of the pres- 
ence of the critical layer (see below). The solution for the 
wavespeed c is assumed of the form 



c — co + eci + ■ 



(41) 



tains only the azimuthal velocity contribution because the 



® With respect to the other scalings we have called on in this 
work, we shall formally require that e <^ 5 <^ e <^ 1. 

Note, however, that we show in Section 4.3 that the results 
obtained in this section, including the calculation of the growth 
rates, would be essentially unchanged had we imposed, instead, 
that the radial gradient of the potential vorticity be zero at the 
starside boundary. 
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and a similar series of the form 
n = Ho + efii + • ■ • 



(42) 



We know apostcriori that the critical layer will generate a 
solution proportional to e In e and this is why such a term 
appears in the above expression. Thus at 0(1) we find 
that 



must be controlled by considering a boundary layer calcu- 
lation in and around the critical layer. Note that since the 
vertical vorticity is proportional to 9^111, the critical layer 
will appear as a vortex sheet. The formal presentation of 
this solution, including the region of validity and expres- 
sion of the boundary conditions at this order is presented 
in Sectionll3.il 



0. 



(43) 



To this order the no-flow boundary condition at a; = 
amounts to 



= codxflo + -flo, 

q 



at X = 0. 



(44) 



The solution to this equation which decays as a; ^ oo is 

Ho = Aoe-''^^ (45) 

where Aq is an arbitrary amplitude. Using this solution in 
boundary condition (|44p amounts to selecting cq, which is 



Co = 



k^q qF^(3 



(46) 



This says that the wavespeed is real and positive, which 
means in this case that there will be a critical layer in the 
domain, i.e. &t x = Xc = cq (see below). 
At O (e) we find the equation 



{c^-x){dl-kl)n. 



Mdl~kl)tl,~t{dl~P^~2q)dl\l,. 



(47) 



Using (j45p for Hq and dividing the equation hy cq — x we 
find the more transparent form 



i{x ~ Co) 



Cid{x - Co). 



(48) 



In which we have defined the parameter A = fc^ (fc^ ~ P'^ — 
2q) . The delta function appearing with the arbitrary coef- 
ficient Ci is a formal device used to signal the presence of 
a critical layer in the flow (Case, 1960, and more recently, 
Balmforth & Piccolo, 2001, Balmforth et al., 2001). The 
coefficient will result from the critical layer analysis (see 
below) which follows the procedures found in similar in- 
vestigations (e.g. Stewartson, 1981, Balmforth & Piccolo, 
2001). In practice it means that we must separately de- 
velop solutions to III on either side ofx = Xc- Before doing 
so let us observe the reason why this series expansion fails 
near Xc by developing the solution to Hi in the vicinity of 
this point. An indicial analysis shows that 



til ^ i ^AAqc [{x - Co) ln(x - cq) - {x - Cq)] -f 
5:^111 ^ r iAAoe-2/« ln(x - cq) + ■ ■ ■ 



1 



X - Co 



(49) 



These expressions show, especially that for that the 

solution begins to breakdown (i.e. break order) when the 
quantity |a; — co| starts to approach 0. This divergence 



4.2.2. Critical layer calculation, matching, and growth 
rate and analysis 

As we demonstrated above, the solutions begin to break- 
down in the vicinity of the critical layer which are those 
places where x — Xc begins to get small. We must therefore 
reexamine (1341) in this zone and to this end we define a 
new inner coordinate as 



,1/3^ = 



(50) 



According to this new coordinate (|34p is reexpressed as 



n = £2/3 



2q)dl 



n 



,2/3 



Cl 



91 



n. (51) 



We intorduce a series expansion for the solution to ([51 
by writing 



n = Hn 



/3 



,2/3n 



2/3 



/3 



(52) 



The remainder of this calculation including the matching 
of the inner and outer solutions and the determination of 
the growth rate ci has been relegated to Appendix IB. 21 
We note that the term proportional to e In e is needed for 
matching purposes as the inner solution is extended out of 
the critical layer (for details see the full exposition in the 
Appendix). We find that the correction wavespeed obeys 



Cl 



^e4/9(fc2 _^2_2g) 



1 



-i/q 



Ei 



q 



The growth rate of this mode, i.e. a 
proportional to 



(53) 



-egfclm(ci), is 



cr = Re-M (F/-l)/3^-2g 



l + ^,-4/,Eif4 



(54) 



where we have restored the definition of e in terms of Re, 
and where we have replaced k^ accordingly with F^p. The 
term inside the square brackets is always greater than zero 
for q > and limits to 3/2 as q 0|j Therefore there is 
growth when F-^ > 1 -I- 2q/P'^. 

The instability emerges from the inner boundary due 
to application of the no-normal flow condition but it will 



Note, however, this limit breaksdown the QHSG approxi- 
mation and is not considered. 
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be at the critical layer where evidence of it appears in 
the form of a pronounced PV-anomaly. The radial width 
of this vortex zone is proportional to e^/^^ the size of the 
box and the amplitude will be t^l'^ times the leading order 
perturbation pressure field (see the end of Appendix IB. 2p . 
Restoring units and recalling that the shearing-box has 
been scaled according to the thermal scale height of the 
disk H , the radial width of this vortex zone, Ai?, is 

AR - a^^^H. 

Further analysis of these solutions together with the 
Reynolds-Orr equation ([55)1 and the definitions (|39ll40p 
shows that to lowest order in Rc^^ 

Stress = Re~^2kFP^F^ + O (Re-') , 

= Re-^kFP^{F^ + 1 + 2q/P^) + O (Rc-2) (55) 

so that 

i^strcss - E^isc = Rc-^kpf^^iFf - 1 - 2q/l3^) + O (Re-2) . 

A general evaluation of the Reynolds stress term £^shcar 
shows that it contributes first at order Re~^ as well, i.e. 

i^shcar = ^(^~^E^lLr + ^ (Re"') , 

(1) _ 



Pressure: n 



PV-anomaly: Q 



shear 



-^shear ^^^^ Order correction to the azimuthal/radial 

velocity correlations. Its details may be worked out but we 
leave it here in general form in order to make the following 
argument. The Reynolds-Orr equation evaluated for this 
set of boundary conditions is to lowest order in Re^^, 

dE „ ^ ^ r, ■ p\F^-l-2q/p^) 



V 



{uiVq + uoVi + c.c.)d 



dt 



kpE, 



(1) 



shear 



2a8 = Re" 

£ = [1^2 + 1] > 0, 
plus a correction which is O (Re^'^) 



1 

(56) 

In the case of instabil- 

£^visc < 0. 



ity, i.e. > 1 + 2q/f3^, it follows that Estr. 
Given ([M]) taken together with ([55)1 we find that under 
conditions of linear instability it also follows that 

P\f! - 1 - 2q/P^) X 



kpE^l, 



1 



1 



-f: 



,-4/9 



1 



(.57) 



4 V V 9 \q. 
Let us reflect upon this for a moment: with these boundary 
conditions it is always the case that i?stross — ^'viscous > 
when F^ > 1 + 2q/P'^. The combined action of the star- 
side viscous stresses and the domain integrated viscous 
lossess still results in a net transfer of energy into the 
domain. Because the term within the square brackets in 
([57]) is always positive, these conditions promote the type 
of corrections (at order Re~^) to the azimuthal and radial 
velocity profiles such that a positive global correlation be- 
tween them emerges and, hence, resulting in E^l^^^^. > 
when F^ > l + 2q//3'^. We conclude that with these bound- 
ary conditions the instability observed is fed both by the 
energy entering the domain from the starside boundary as 
well as by the energy extracted by the shear due to the 
resulting order Re~^ velocity profiles. 




(d) Stress-fre# 
at X = 
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Fig. 1. A comparison of eigenfunctions for a variety of 
starside boundary conditions. In all plots f3 — 2, Fc — 1.2 
and e = 10^^. The pressure eigenfunctions II are shown on 
the left panel of plots while the potential vorticity eigen- 
functions are shown on the right panel. The predicted 
growth rates are also quoted, (a) Zero potential vorticity, 
(b) zero potential vorticity gradient, (c) rigid boundary, 
(d) stress-free and (e) zero-pressure fluctuation. . 
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Fig. 2. Like Fig. [T] except e = 10' 
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4.3. Finite radial domain investigations 

In this section we consider the normal-mode solutions of 
occurring on a finite radial domain where Xout = 1- 
All solutions are computed numerically using a Newton- 
Raphson scheme on a Chebyshev grid of anywhere from 
33 to 129 points - higher resolution is needed for smaller 
values of e. All numerically generated solutions are nor- 
malized so that fldx = 1. Because this is a fourth order 
system we must specify four boundary conditions. In all of 
the following calculations two of the boundary conditions 
will be that 

n = 0, Q = 0, at X = 1, (58) 

in other words, that the pressure fluctuation and PV- 
anomaly are zero on the farside boundary. According 
to the normal-mode PV-anomaly, Q = d^Il — kpYl, the 
fixed pressure condition implies that at the far boundary 
din = 0. It therefore follows from ^ that the perturbed 
stress expression, dxV, is zero there as well. A third con- 
dition will be that there is no normal-flow on the starside 
boundary (as before), i.e. 

u = 0, at x = 0. (59) 

For the remaining starside boundary we shall explore four 
different conditions enumerated in the corresponding sub- 
sections below. The most physically plausible viscous star- 
side condition is to set to zero either the azimuthal velocity 
fluctuation or the azimuthal stress fluctuation. We have 
also considered zero PV-anomaly and zero PV-anomaly 
gradient conditions. Although these conditions are less 
physically realistic, they are simpler to interpret in terms 
of the energy arguments developed in previous sections. 

4.3.1. Rigid and stress-free starside boundary 

These conditions translate to requiring v — (rigid) or 
dxV = (stress-free) at 2: = 0. Note that the usage of 
the term "stress-free" to describe the boundary condition 
really refers to the perturbed part of the azimuthal flow 
as being stress-free. Inspection of Fig. [3] shows that there 
are regions in the e-i^e parameter plane in which there is 
instability when rigid or stress-free conditions are imposed 
on the starside boundary. Although the parameter range 
for instability is not as expansive as it is for the other less- 
realistic boundary conditions explored (see both the pre- 
vious section and below) , the results suggest that a rough 
criterion for linear instability is > 1. Furthermore, un- 
der these starside conditions we find that the Reynolds- 
Orr equation (|38p simplifies to 

— r- = -Eshcar — -Evisc, (60) 

at 

since the zero-stress and rigid conditions implies that 
^'stress = 0. Thus wc find that the boundary conditions 
are such that instability is not directly driven by the in- 
jection of energy into the domain due to the (perturbed) 



body stresses. Instead we interpret the instability as a con- 
sequence of the velocity profiles set up by the conditions. 
In other words, since i?visc is positive definite, this pro- 
cess experiences growth entirely due to the extraction of 
energy from the shear as embodied by the domain integral 
term -Eshoar- The character of the eigenfunctions are seen 
by inspecting the corresponding profiles for small values of 
e in Figure [1] In both the rigid and stress-free cases there 
are prominent boundary layers appearing on the starside 
for both the pressure and potential vorticity. The critical 
layer in the potential vorticity also appears here (see the 
inset in Figure [TJ; and [TJl) but is dwarfed by the starside 
boundary layer. 

4.3.2. Zero PV-anomaly 

This boundary condition may be envisioned as the starside 
boundary counteracting any tendency for the development 
of any PV-anomaly there. Although this is somewhat arti- 
ficial, we present here the results of this investigation be- 
cause these boundary conditions give solutions that closely 
resemble those obtained for the calculation on the semi- 
infinite domain calculation. As in the semi-infinite domain 
calculation, instability occurs when > 1 + 2q/P^, and 
it scales as e for e <C 1. Inspection of the eigenfunctions in 
Figure [T] for e ^ 1, especially the profiles for Q, shows that 
(i) the perturbed potential vorticity is strongly localized 
in the critical layers occurring where the real wavespeed 
approximately equals the background flow speed, (ii) an 
additional boundary layer appears at the starside bound- 
ary scaling like e and, (iii) the vorticity in the critical layer 
follows the e-^^^ scaling determined in the semi-infinite do- 
main calculation. We depict in Fig. 2^ a contour plot of 
growth rates as a function of both the Reynolds number 
and the inverse of the Froude number, F^^ for fixed values 
of (3 and k. The vertical axis may be understood as repre- 
senting a positive increase in the wave's speed (see Eq.l46|). 
For the parameters depicted in Fig. (/3 = 2, g = 3/2), 
instability sets in for F^ > •\/7/4. 

4.3.3. Zero radial PV-anomaly gradient 

Requiring no radial gradient of the PV-anomaly on the 
starside boundary is arguably the least physically real- 
istic but we include it here, as in the previous section, 
because it best reproduces the asymptotic result of the 
semi-infinite domain calculation. Like in the previous sec- 
tion, where Q is set to zero there, there is instability when 
F^ > 1, a, critical layer emerges which also scales as e^/'^ 
for small values of e. However, the boundary layer ap- 
pearing near the starside boundary for the calculation of 
Section 14.3.21 vanishes here. Finally, for small values of e 
(i.e. Re^^) the Reynolds Orr expression for these distur- 
bances takes on the same leading form as (|56p in Section 
14.21 and this includes the character of Sgtross- We note 
that growth rates here are nearly identical to the growth 
rates determined in Section 14.3.21 Figure Hb shows the 
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(a) Stress-Free Boundaries 




(b) Rigid Boundaries 




12 3 4 

log(Re) 



Fig. 3. Contours on the F^-Ke plane of Im(— c) for q = 
3/2, k — 1 and (3 — 2. Given the value of fc, according 
to its definition e is written in terms of Re = 3 /{2a): 
(a) Stressfree boundary conditions at x = 0, (b) Rigid 
boundary conditions at a; = 0. Shaded regions indicate 
growing modes. Note that according to (1^5]) the vertical 
axis indicates increasing wavespeed. 

landscape of instability (for f3 ~ 2,q — 3/2) and we see 
that instability also sets in when > 1 + 2q/ 0^ (i.e. here 
for > a/7/4) but that it is bounded above by a more 
complicated function of Re. We note that the critical layer 
becomes harder to distinguish as viscosity (that is, Re~^ 
or e) is made larger. 

5. Discussion and Reflections 

5.1. On the Tollmien Schlichting wave analogy 

T-S waves appear in flows that are are wall-bounded at 
least on one boundary. In the classic analysis done for 
Blasius boundary layers (e.g. Schlicting, 1968) the insta- 
bility is a solution of the 2D Orr-Sommerfeld equation. 
The mechanics leading to instability is understood to arise 



(a) Zero PV-anomaly at x=0 




(b) Zero PV-anomaly gradient at x=0 
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Fig. 4. Same as Fig. [3] except: (a) Zero PV-anomaly at 
a; = 0, (b) Zero PV-anomaly gradient at x = 0. 

from the action of a purely viscous mode interacting with 
a nearly inviscid mode (Baines, Majumdar & Mitsudera, 
1996). The global velocity profiles set up are such that the 
relative phase between the horizontal and (plate) normal- 
velocities promotes extraction of energy from the shear 
which then leads to energetic growth (i.e. Fshcar > 0). 
The 3D instability studied here shares some major simi- 
larities to classic T-S waves: (i) the equation governing the 
dynamics of the potential vorticity modes (j30p has similar 
structure to the 2D Orr-Sommerfeld equation, (ii) the in- 
stability emerges for both no-slip and free-slip boundary 
conditions (but not limited to these). What stands out in 
our minds is that, although classical unstable T-S waves 
come about in a wide variety of background flows exclud- 
ing plane-Couette flow (Baines, Majumdar & Mitsudera, 
1996, Schmid & Henningson, 2000), the instability here is 
present for a plane-Couette type of flow profile. 

The results here compare qualitatively to the results 
of rotating Blasius boundary layers. For example, Yecko 
& Rossi (2004) show that three dimensional modal in- 
stability preferentially emerges in anticylonic rotating 
Blasius flow when the vertical wavenumber of the distur- 
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bance is large in comparison to its streamwise (azimutlial) 
wavenumber (e.g. see Figure 8b in Yecko & Rossi, 2004). 
By comparison, the asymptotic scalings we have imple- 
mented spotlights dynamics characterized by these same 
spatial scale disparities. Thus despite the differences in the 
problems investigated between these two studies (i.e. the 
inclusion of gravity and entropy gradients and the differing 
base velocity profiles) the similarities in the circumstances 
for instability onset suggest that such processes may be 
more general in environments like disks - especially near 
the disk-star boundary. 

5.2. On the assumed constancy of g and N 

To make the analysis we have exposed here analytically 
possible we assumed that the vertical component of grav- 
ity and the Brunt- Vaisala to be constant with respect 
to the disk vertical coordinate, ([51]) which permits us to 
assume separable normal-mode solutions. On the other 
hand, real disks (including their small sections) are char- 
acterized by vertically varying values of g and N and this 
means that, in general, one cannot assume separable so- 
lutions in z and x, in particular. 1^1 We have checked that 
the results obtained in the limit where the viscosity pa- 
rameter is small (i.e. e ^ 1) still holds when g and N are 
taken to be correctly z dependent. Accordingly we have 
repeated the asymptotic calculation described in Section 
4.2 where, in addition, we restricted our attention to finite 
vertical domains by imposing either velocity or pressure 
conditions on the vertical boundaries. Such disturbances 
will be characterized by vertical overtones labeled by an 
overtone wavelength /3„ - which should be thought of as 
being analagous to the vertical wavenumber (3 we assumed 
in Section 4. The asymptotic calculation shows that insta- 
bility sets in so long as a vertically weighted Froude num- 
ber, Ff{Pn), approximately exceeds one. The calculation 
is far more lengthly and does not add any new qualitative 
details to the one presented in this work and it is for this 
reason we have omitted it from the current exposition and 
we will expand upon it in a future study. 

5.3. Relationship to the viscous pulsational instability 

In an axisymmetric study of an shearing sheet section 
of an accretion disk of constant temperature, Latter & 
Ogilvie (2006) argue that the VPI (Kato, 1978) is most 
likely to manifest itself through the destabilization of an 
even structured f-mode. The disturbances become unsta- 
ble because the viscous perturbations transfer energy from 
the shear into the acoustic mode through the perturbed 
stress r^y. Because f- modes are characteristically inertial- 
acoustic waves, they are the likely candidates for this in- 



In other words by assuming g and A'^ to be constant we are 
able to assume solution form 11 = X{x)Z{z)T{t)Y{y) while if 
g and A*' are z-dependent one can (at best) assume a solution 
in the form 11 = Tl{x,z)T{t)Y{y) where n(a;,2:) is the non- 
separable structure function. 



stability since their vertical structures are the simplest 
which, in turn, result in minimizing dissipative losses. We 
observe that T^y is proportional to the pressure fluctua- 
tion. This fluctuating stress affects the evolution of the 
horizontal velocity perturbation by extracting/adding en- 
ergy into the disturbance. As the horizontal velocity dis- 
turbances are not in general in phase with the pressure 
fluctuations, especially for inertial/acoustic modes, the 
possibility for overstability is manifest. However, the PV- 
anomalies examined here are distinct from inertial-gravity 
and inertial-acoustic modes (e.g. Ogilvie, 1998, Tevzadze 
et al., 2004). The horizontal velocity perturbations of PV 
anomalies are proportional to the radial gradient of the 
perturbation pressure. Since the latter quantity is propor- 
tional to a decaying exponential (i.e. ^ e~^f''^), it follows 
that the phase between T!j.y and the energy in the PV- 
anomaly will be tt radians out of phase with each other 
(see Section 4.1). Thus the fluctuating stress T^j, behaves 
to stabilize a PV-anomaly. resulting in a stabilizing rela- 
tionship as we observe in Section 4.1. As such we under- 
stand the destabilization of the PV-anomaly as as being 
distinct from the instability leading to the VPI. 

5.4. Summary and Implications 

For cold disk systems, i.e. those in which magnetic ef- 
fects are not active, the prospects of identifying insta- 
bility mechanisms appear to be far from exhausted. We 
have tried to argue that certain previously considered 
non-magnetic instabilities need not be discarded as can- 
didate mechanisms driving activity for disks. Indeed the 
SRI and PPI instabilities, which emerge as the interaction 
of edgewaves along cylinder/channel walls, could in princi- 
ple operate in real disks so long as there exists, in general^ 
interacting waves propagating separated from each other 
by a wave-evanescent region. We have demonstrated here 
another possible mechanism - that the existence of disk 
analogs of unstable ToUmien-Schlichting waves could also 
manifest themselves in real disk systems. We have car- 
ried out the calculation within a model shearing box in 
which we have imposed a single boundary on one side. 
True disks have boundary layers separating stars from the 
disk which are probably far more complicated (Regev & 
Bertout, 1995) than the model we have presented here. 
Nonetheless, far from being a proof, we have demonstrated 
in this asymptotically simplified model that such a dy- 
namical processes is, at least, feasible. It is no stretch 
of the imagination to suppose that analogous unstable 
waves may exist near the vicinity of the star-disk bound- 
ary layer. We add a final refiection. Classical unstable T-S 
waves emerge in flows with compliant boundaries showing 
that such instabilities are robust and persist even if the 
boundaries have a certain amount of elasticity to them 
(Carpenter & Garrad, 1985) - although compliant walls 
delay the onset of instability to higher Reynolds num- 
bers. As a star-disk boundary is probably not a rigid body 
transition, it would be beneficial to investigate and/or 
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model these disk analog T-S waves by considering starside 
boundary conditions that are appropriately compliant as 
well. 

The perspective we have adopted therefore can be bro- 
ken down into two parts. The first is that (as in the recent 
studies of Kleiber & Glatzel 1999, and Latter & Ogilvie 
2006) we have expanded the exploration of the possible 
destabilizing role viscosity can play. Viscosity does not al- 
ways have to stabilize disturbances as there are velocity 
profiles, dictated by boundary conditions, wherein desta- 
bilization occurs counter to one's usual physical intuition. 
The second is that experimentation with boundary condi- 
tions, even within the context of the shearing box environ- 
ment, followed by concerted effort toward understanding 
and clarifying their effects is a worthwhile endeavor given 
our lack of complete knowledge about the boundaries of 
real disk systems (a situation which is strongly contrasted 
by what is encountered in laboratory/terrestrial flows). 

If turbulent stresses in cold disks are driven by the MRI 
resulting in effective values of a 10^'^ — 10^'* (Ogilvie, 
2003, King, Pringle & Livio, 2007), then T-S waves hke 
the sort here could emerge as a secondary instability. This 
possibility is made manifest because the T-S wave insta- 
bility grows faster in proportion to the value of a until 
about a value of 0.1. 

6. Acknowledgements 

The authors are indebted to the valuable comments and 
suggestions of the anonymous referee. The authors would 
like to thank the Israeli Science Foundation for mak- 
ing this research possible. OMU also acknowledges that 
this research was partly supported by BSF grant 2004087 
and ISF grant 1084/06. OMU also thanks the Dead Sea 
Regional Council and the Ein Gedi Kibbutz for their hos- 
pitality and Phil Yecko for suggestive conversations. 

References 

Afshordi, N., Mukhopadhyay, B., & Narayan, R., 2005, ApJ, 
629, 373 

Baines, P. G., Majumdar, S., & Mitsudera, H., 1996, J. Fluid 
Mech. 312, 107 

Baines, P. G. , & Mitsudera, H., 1994, J. Fluid Mech. 276, 327 
Balmforth, N.J., J. Fluid Mech., 387, 97 
Balmforth, N.J. & Piccolo, C, 2001, J. Fluid Mech., 449, 85 
Bahnforth, N.J., Piccolo, C, & Umurhan, O.M., 2001, J. Fluid 

Mech., 449, 115 
Balbus, S.A., 2003, ARAA, 41, 555 
Barranco, J.A., & Marcus, P.S., 2005, ApJ, 623, 1157 
Bayly, B. J., Orszag, S. A., & Herbert, T., 1988, Ann. Rev. 

Fluid Mech., 20, 359 
Bodo, G., Tevzadze, A., Chagelishvili, G., Mignone, A., Rossi, 

P., & Ferrari, A. 2007, A&A, 475, 51 (Bodo et al. A) 
Bodo, G., Chagelishvili, G., Murante, G., Tevzadze, A., Rossi, 

P., & Ferrari, A. 2007, larXiv : 0705 ■ 3474t f 1 Fastro : phi 

(Bodo et al. B) 

Carpenter, P. W., & Garrad, A. D. 1985, J. Fluid Mech., 155, 
465 

Case, K.M., 1960, Phys. of Fluids, 3, 143 



ChagelishviU, G.D., Zahn, J.-P., Tevzadze, A. G., & 

Lominadze, J.G., 2003, A&A, 402, 401 
Charney, J. G. & Stern, M. E. 1962, J. Atmos. Sci., 19,159 
Davies, H.C., & Bishop, C. H. 1994 J. Atm. Sci., 51, 1930 
Doering, C.R, Spiegel, E.A., & Worthing, R.A., 2000, Phy. 

Fluids, 12, 1955 
Drazin, P.G. & Reid, W.H., 1984, Hydrodynamic Stability, 

Cambridge. 

Dubrulle B., Marie L. , Normand Ch., Richard D. , Hersant 

F.& Zahn J.-P. 2004, A&A, 429, 1 
Frank J., King A.R. & Raine D.J. 2002, Accretion Power in 

Astrophysics. Cambridge Univ. Press, Cambridge 
Goldreich, P., Lynden-Bell, D., 1965, MNRAS, 130, 125 
Goldreich, P., Goodman, J., & Narayan, R. 1986, MNRAS, 

221, 339 

Hayashi, Y.-Y. & Young, W.R. 1987 J. Fluid Mech., 184, 477 
Hoskins, B. J., Mclntyre, M. E., & Robertson, A. W. 1985 

Quart. J. Roy. Meteor. Soc, 111, 877 
loannaou, P.J. & Kakouris, A., 2001, ApJ, 550, 931 
Ji, H., Burin, M., Schartman, E., & Goodman, J. 2006 Nature, 

444, 343 

Johnson, B.M. & Gammie, C.F., 2005 ApJ, 635, 149 
Joseph, D.D., 2003, J. Fluid. Mech. ,479, 191 
Kato, S., 1978, MNRAS, 185, 629 

Kersale E., Hughes, D. W., Ogilvie, G. I., Tobias, S. M., & 

Weiss, N. O. 2004, ApJ, 602, 892 
Klahr, H.H. & Bodenheimer, P., 2003, ApJ, 582, 869 
Kleiber, R. & Glatzel, W., 1999, MNRAS, 303, 107 
King, A. R., Pringle, J. E., & Livio, M. 2007, MNRS, 376, 1740 

(KPL07) 

Kluzniak W. & Kita, D. 2000, Three-dimensional structure of 

an alpha accretion disk, astro-ph/0006266 
Latter, L. N., Ogilvie, G. I., MNRAS, 372, 1829 
Lesur, G. & Longaretti, P-Y., 2005, A&A, 444, 25 
Li, H., Finn, J. M., Lovelace, R. V. E. & Colgate, S. A, 2000, 

ApJ, 533, 1023 

Lions, P.L. 1993, Limits incompressible et acoustique pour des 
fluides visqueux, compressible et isentropique, C.R. Acad. 
Sci. Paris Ser. I Math, 317,1197. 
Lithwick, Y. 2007 ApJ, 670, 789 (Lithwick 2007 A) 
Mukhopadhyay, B., Afshordi, N., & Narayan, R., 2004, ApJ, 
629, 383 

Ogilvie, G. I., 1998, MNRAS, 297, 291 
Ogilvie, G. I., 2001, MNRAS, 325, 231 

Ogilvie, G. I., Proctor, M. R. E., 2003, J. Fluid Mech., 476, 
389 

Papaloizou, J.C.B., & Pringle, J.E. 1984, MNRAS, 208, 721 
Petersen, M. R., Julien, K., & Stewart, G. R. 2007 ApJ, 658, 

1236 (Petersen et al. A) 
Petersen, M. R., Stewart, G. R., & Julien, K. 2007 ApJ, 658, 

1252 (Petersen et al. B) 
Regev, O. & Bertout, C. 1995, MNRAS, 272, 71-79 
Richard, D., 2001, These de doctorat, Universite Paris 7. 
Richard, D., Zahn, J.-P., 1999, A& A, 347, 734 
Schlichting, H. & Gersten, K. Boundary Layer Theory, 2001, 

Springer, p 94. 
Sakai, S., 1989, J. Fluid Mech. 202, 149 

Schmid, P. J. & Henningson, D.S., 2000, Stability and 

Transition in Shear Flows, Springer 
Shalybkov D. & Riidiger G., 2005, 2005, A&A, 438, 411 
Spiegel, E. A. & Veronis G., 1960, ApJ, 131, 442 
Stewartson, K., 1981, IMA J. Appl. Math., 27, 133 
Tevzadze, A.G., Chagelishvili, G.D., Zahn, J.-P., Chanishvili, 

R.G., & Lominadze, J.G., 2003, A&A, 407, 779 



14 



Umurhan, O.M. 2006, MNRAS, 365, 85 

Umurhan, O.M., Nemirovsky, A., Regev, O., & Shaviv, G., 

2006, A&A, 446, 1 
Umurhan, O.M & Regev, O. 2004, A&A, 427, 855 
Yecko, P.A. 2004, A&A, 425, 385 
Yecko, P.A. & Rossi, M. 2004, Phys. Fluids, 16, 2322 

Appendix A: Scaling arguments leading to the 
viscous large shearing box equations 

The derivation of the VLSB equations follows the proce- 
dure executed in Umurhan & Regev (2004). The dimen- 
sional equations of motion in cylindrical coordinates in a 
frame of constant rotation are, in component form, given 
by the following 

~o ~ \ BP 

^-nlr-2noU,-^ = --^-^ + N, (A.l) 
at r p Br Br 
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pN^ 



1 Brt, 



1 Brt 



r Br 

1 BrH^ 

2 Qj. 

1 drt_ 
r 



<t"t> ^ 
r 

dt. 



Bz 



r 



r Bcf) 



Bz 



Bt 



Br r Bcj) Bz 
along with the viscous stresses t. 

t = 



du, ~ u.u 
at r 

~di 



1 BP la* 

pr Bcj) r B(j) 
1 BP 9* 



t = < _ = 




2,(1?^ + ^ 



p Bz Bz 

the equations of mass continuity and entropy 



dp ^ /I Bru^ ^ 1 Bu^ ^ Bu^ 
dt ^ \ r Br r B(b Bz 



0, 



with the operator definition 



B 



d B 



dt Bt '' Br r 



' Bz- 



(A.4) 
(A.5) 



t^^ 

in which 

Vr • U = 



r 
Bu 
Bz 



2 

Bz 



r 

iBu 

r B(j) 



+ c 



3^ 



Vr ■ U, 



(A.7) 
(A.8) 
(A.9) 

(A.IO) 

(A.ll) 

(A.12) 
u,(A.13) 
(A.14) 
(A.15) 



1 Bru 1 Bu, 



Br 



Bru^ 
Bz 



u^,u^, are the radial (r), azimuthal (0) and vertical (z) 
velocities as observed in the rotating frame. The rotation 
rate fio is set to the rotation at some fiducial radius r = 
Rg . The curious notation on ii^ is meant to indicate that 
the azimuthal velocity observed in the laboratory frame 
(denoted by u^) would be related to its velocity in the 
rotating frame by = il^r ~\- u^. The entropy is defined 
to be S = h\Pp~^ in which is the specific heat at 
constant volume, 7 is the ratio of specific heats, that is 
7 = Cy/Cp in which Cp is the specific heat at constant 
pressure. Q is a heat function representing some sort of 
non-adiabatic processes that are relevant to the disk. 

The gravitational potential 5* is written in the follow- 
ing unusual form in order to effect some generality when it 
comes to the shear it induces upon the steady state flow. 



Because the turbulent viscosity within a disk is presumed 
to be driven by some sort of shear process (either the 
MRI or a subcritical hydrodynamic transition) the bulk 
viscosity C, will be set to zero hereinafter. According to the 
classic proposal of Shakura & Sunyaev (1973), the shear 
viscosity is parametrized as 



2P 



a 



(A.16) 



where Q, = f2o(-Ro/'')'^ is the steady-rotation rate induced 
by the generalized potential (|A.6|) . When we consider 
Keplerian flows we write Vl^ = il{q — 3/2) (and see be- 
low). The parameter a is a tunable order 1 quantity. 

We proceed from here onto non-dimensionalization. 
We shall consider dynamics as taking place in a small 



region around the point r = Rot4' = 



= 0, which 



we shall refer to as the box. We let e measure the non- 
dimensional size of this box, i.e., 
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For the realistic Keplerian flow profile q = 3/2 and 
fl^ = GM^./R'^, where is the mass of the central ob- 
ject. However, r2„ is in general a function of q and only 
when q = 3/2 it is to be interpreted as the "Keplerian" 
rotation rate. 



Ro^<z<R,-. 



This motivates us to scale all the directions by R^ and to 
define the nondimensional coordinates 



ex = 
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where it is understood that x, y, z are now taken to be 
order 1 quantities. 

Furthermore we suppose that all density quantities are 
scaled by the reference density p, pressure quantities are 
supposed similarly scaled by pc^ where Cg is the dimen- 
sional scale of the sound speed of the box. In parallel with 
this speed is the local rotation speed of the box around 
the central object, = i^gRg- when q — 3/2 this speed is 
sometimes referred to as the local Keplerian speed of the 
disk. The fundamental ansatz of cold thin disk theory is 
that the ratio of these two quantities is small. We, in fact, 
identify 



V ' 

* n 



(A.18) 



which is the classic parameter measuring the thinness or 
"coldness" of the disk (Shakura & Sunyaev, 1973 and see 
recently Umurhan et al., 2006, for recent exploitations of 
this parameter). This, in turn, is a measure of the disk's 
vertical scaleheight H = eRg. We also note here that by 
equating the size of the box to the " coldness" of the disk 
means we are here looking at the viscous analog of the so- 
called "Large Shearing Box" equations formally developed 
in Umurhan & Regev (2004). 

Thus we propose that all velocities observed in this 
rotating frame are scaled by the soundspeed, meaning to 
say that 



where u',w',v' are order 1 nondimensionalizations of the 
corresponding velocities. Note that in these scalings we are 
saying that all velocities observed in the rotating frame are 
order e smaller than the basic rotation speed of the box Vg . 
This seemingly obvious point is emphasized because the 
magnitude scale of the accretion and meridional flow in- 
duced by the turbulent viscosity is order smaller than 
Vg (Shakura & Sunyaev, 1973, Kluzniak & Kita, 1999). It 
means that to the order to which matters are considered 
here, that is as far as the generalized "box" formalism is 
concerned (see below), the effects of accretion and merid- 
ional flow are absent. . 

Time and all advective derivatives are scaled accord- 
ing to the local rotation time of the disk, i.e. Clg. Before 
putting these scalings into the governing equations we note 
that 

^l^-^-'^^lRo [eqx + 0{e^)]. 

Taking into account all of the nondimensionalizations, 
along with writing p pp' ^ P pc^p' and understanding 
that p',p' are the nondimensionalized density and pressure 
quantities, we find that the equations of motion (|A.HIA.3|) 
become 
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N'^, N'^, N'_^ are as they appear in the text. To reiterate, 
the effects grouped in the O (e) terms consist of curvature 
effects, higher order corrections due to the central poten- 
tial, and turbulent viscosity induced accretion/meridional 
flow. The equation of continuity and entropy conservation 
also become 

dp' , ,dp' ^,dp' ,dp' 
dt dx dy dz 

To obtain the viscous large shearing box equations pre- 
sented in the text, i.e. ([T][5]), we make the following iden- 
tifications and assumptions 



— Drop all terms O (e) and higher from (|A.19IIA.23|) . 

— Everywhere write w' — —qx + w' in (|A.19IIA.23)l in or- 
der to eliminate the sole x term on the RHS of (jA.19[) . 
The expression —qx is the background shear and will 
be felt by perturbations. 

— Write the density and pressures as being comprised 
of a steady portion and a time dependent perturbed 
portion, i.e. 

P' Pt{x,z) + p{t,x,y,z), 

— Specifically flag all Coriolis-like related effects with the 
symbol 17^ = 1. 

Appendix B: Semi-infinite domain calculation 
details 

B.l. Outer solution completion 

As a result of this breakdown we write the solutions to 
this order according to which side of Xc one is on. Formally 
then we say 

Hi = n^"' = A^^^e-'^^-^ + Sj^^e'^^-^ + (B.l) 



for < a; < — e_, and 

iii = n(+)=4+)e-'^-.--Kni+\ 



(B.2) 



for X > Xc + S-f-. The particular solutions n^p^ satisfy 
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with solutions given by 
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where Ei(a;) is the Exponential Integral function 
(Abramowitz & Stegun, 1972). We define 



^ip' (a^c ) = lim n^p 
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An analysis of these solution forms for n^p'' shows that 
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In formally writing these solutions we purposely avoid the 
region described by 

Xc — £_ < X < Xc + £+, 

which we will refer to as the critical layer. The size of this 
bounding region is such that 0(1) > 0(e±) > 0(e). 
It will be shown that the solutions on either side of this 
region will asymptotically match at lowest non-trivial or- 
der to the solution emerging from the critical layer itself. 
Finally, the u = boundary condition at 2; = may be 
read from (|36|) after making use of the above relationships 
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B.2. Critical layer calculation and matching 
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The calculation will be facilitated if we consider the evolu- 
tion of the residual potential vorticity quantity n through 
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The polynomial terms in the above expression are the first 
five terms of the series expansion of the lowest order outer 
solution, Ho expressed in the vicinity of the critical layer. 
We may now rewrite ([?T|) instead in terms of II, 
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where A is as defined in the text. At 0(1) we have simply 

1 



Cd^ni = 0, with C = - . (B.8) 

Homogeneous solutions of the operator Cd^ of these in- 
volve integrals of Airy Functions associated. However, 
because these are functions of complex arguments these 
two solutions are rejected because they show exponential 
growth as ^ ^ ±00 (e.g. Bender & Orszag, 1978). All 
homogeneous solutions involving these Airy functions are 
rejected henceforth. Thus we assume the following solu- 
tion expansion for H 

H = eHi + e^/^n4/3 + • • • 

The order e solution will be Hi — Ai. There is also a so- 
lution at this order proportional to ^ but it is rejected on 
account of the fact that it will not match any correspond- 
ing solution from outside. The equation at order e^^^ is 



CUi/3 = -jAAo. 

with Aq = Aoe"2/«. With solution 

114/3 = ^4/3'? + 114/3 

In which 9|ll4/3 = S and where 
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Operating on S by £ followed by an integration by parts 
verifies that II4/3 is the particular solution of (|B.9p . We 
note immediately that S is bounded for all real values of 
X by A^qF (I) since 



/ exp (iw^ - icj^)dtj < / exp (-icj^)(iw = F I 
Jo Jo 

It is instructive to note the expression 
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where A± > will be related to the expressions of e± asso- 
ciated with the outer solutions (see next section). This re- 
lationship plays the role connecting the two homogeneous 
outer solutions. When A± » an asymptotic evaluation 
of the integral shows that 



/ Ed^ - -Aoe-^/'' [In - In A_ + in] 

J -A * 
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We note that there exists a phase factor proportional to 
—in that arises from this operation. This shall be explic- 
itly referred to in the next section. One final note: in the 
process of matching to the outer solutions the expression 
e^/'^Il4/3 produces terms proportional to elne in the outer 
region. In order to have this appropriately matched (in 
this case, cancelled) is the reason why there is an O (e In e) 
term in the critical layer expansion (|52p . 
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We note that the PV-anomaly arising from this critical 
layer appears at order e^/'^ when viewed in the unstretched 
coordinate frame. This is because the first non-trivial con- 
tribution arising to the potential vorticity perturbation 
from this zone is 

Q = (5^-4)(no + .--) 

^ .4/^9^114/3 + •■• 

Thus, while the size of the critical layer zone is order e^/"^ 
the magnitude of the PV-anomaly will be order e^/^ the 
leading order pressure perturbation. 

B.3. Matching 

We must reexpress the outer solution in a "small" vicinity 
of the critical point Xc- We consider first the solutions 
approaching from below, that is x ^ x~ . We have 
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while when approaching this point from above, that is as 
X — > we find 



1 kp{x — Xc) ~^ h^t-ix ~ Xc)^ 
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Now we do the same to the critical layer solution. We 
restore the inner coordinate ^ in terms of a; — Xc and take 
the limit of small e (as is standard practice in boundary 
layer theory, Bender & Orszag, 1999) revealing 

n = io + kfBj^/3{x - Xc) + ^Aaklix - Xc)^ 
+ iSi/3fc^ {x - Xc)' + j^Aokl {x - Xc)^ 
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Identifications are made respecting powers of e and x — Xc- 
The 0(1) matchings are straightforward since the solu- 



tions coming in from the left and from the right of Xc are 
the same: 



An 
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while the remaining two terms, proportional to [x — Xc)' 
and {x — Xc)'^ are satisfied given the above assignments in 
(|RT7l) . At O (e) we find first that 

e'V,A[~^ + e^/^s'-) + iv[-\xc) = ii = 

because n^p'' (a;~) — Hi^'' (a;^) the above relationship im- 
plies 



(B.18) 



To complete the matching we must prepare the final term 
e"'/^Il4/3. It is asymptotically correct to do this by con- 
sidering the derivative of these terms up to the bounding 
region of of the critical layer, i.e. 1 ^ |a; — Xc| ^ 0. Thus 
we consider the derivatives as one approaches this zone 
(and measured by the coordinates e±) as 
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that is, approaching Xc — £_ from below and 



(B.19) 



eA 



A 



— 



,2/9 



lne+-efc^Aj+^e'2/g^(B 20) 



- a:— »a:c-t-£+ 

that is, approaching Xc + from above. Subtracting the 
two expressions gives 
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Now we must match this to the corresponding expression 
emerging from the interior of the domain. In other words 
we require similarly that 
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in which (i) the transition in orders of e occurs because 
of the change of variables from x to ^, (ii) we have iden- 
tified A± <-> e±/e^/^ and (iii) used (jB.12[) in writing the 
above expression. We noted that the leading behavior of 
the integral for e^/e^/^ ^ 1 (i.e. for the matching zone) 
is 
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Equating the RHS of (|B.2ip and (|B.22|) and making use 
of the asymptotic form (jB.23[) we see that the offending 
logarithmic terms cancel leaving, 
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AQ--m = 74oA7r.(B.24) 
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The meaning of this relationship is clear - the presence 
of the singular layer in the flow means that the homo- 
geneous outer region solutions must show a jump in their 
derivatives (in proportion to the RHS of the above expres- 
sion). Another way to interpret this is to recognize that 
this jump corresponds to the presence of a vortex sheet at 

The complex wavespeed may be now obtained from 
l|B.18|) . (|B.24p and the boundary condition (|B.6p . We note 
that for there to be a non-trivial solution to (jB.lSp and 
(jB.24p the following relationship, 
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must be satisfied. This is the solvability condition. The 
second matter we note is that 
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Thus, the solvability condition (|B.25p together with the 
above expressions used in (jB.6p combine to give 
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The form quoted in the text follows from restoring the 
definition of A into the above expression. 



